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DIGEST 


Flashing  of  agents  during  the  dissemination  process  *  educes  the 
effectiveness  of  chemical  munitions.  Aiterou rul~  'v#  detonation  pro¬ 
ducts  of  the  explosive  charge  appears  to  trigger  flashing,  i.e.,  the 
ignition  of  agent-air  mixtures.  The  chemical  composition  of  the  explo¬ 
sive  charge  detonation  products  controls  the  lgnltlblllty  of  product/alr 
mixtures.  If  product/air  mixtures  cannot  burn,  one  of  the  main  causes 
of  flashing  will  have  been  eliminated. 

At  present  detonation  product  compos 1 tiot. i  cannot  be  obtained  ex¬ 
perimentally.  They  can  only  be  estimated  by  aeoretlcal  calculations , 
some  of  which  are  based  on  Inadequate  or  even  invalid  assumptions.  Con¬ 
sequently,  the  objective  of  this  effort  was  to  choose  the  most  plausible 
of  the  currently  available  theoretical  calculations  and  to  use  it  in  es¬ 
timating  product  compositions  of  likely  dissemination  charges.  Addition¬ 
ally,  such  calculat4ons  provide  estimates  of  other  dotonation  parameters, 
such  as  temperature  and  pressure, which  rav  also  be  important  in  some 
facets  of  the  dissemination  process 

The  computational  scheme  adopted  was  RUBY.  It.  was  used  to  estimate 
detonation  parameters  for  Comp  B,  Comp  B/oxaaide,  amatols,  HMX,  RMX/AP, 
and  HMX/AN.  Based  on  these  calculations  explosite  charges  for  diagnostic 
tests  of  agent  degradation  are  suggested. 
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I  INTRODUCTION 


A.  Objective  and  Outline  of  Approach 

Chemical  agents  may  degrade  during  their  dissemination  by  high 
explosive  munitions.  One  of  the  more  spectacular  forms  of  agent  de¬ 
gradation,  though  not  necessarily  the  only  one,  is  called  "flashing." 

As  the  term  implies,  flashing  refers  to  the  formation  of  a  fireball  of 
burning  agent  and  high  explosive  detonation  products.  Clearly,  flashing 
or  other  degradation  reduces  the  effectiveness  of  a  chemical  munition. 
Consequently,  the  objective  of  SRI 's  overall  p^ogras  is  the  elimination 
of  agent  degradation.  The  particular  phase  of  the  program  reported 
here  deals  with  the  early  stages  of  the  dissemination  process,  namely, 
the  interaction  of  high  explosive  detonation  products  with  the  agent. 
Calculation  of  detonation  parameters  such  as  product  composition  and 
temperature  could  lead  to  (1)  a  rational  understanding  of  the  role  of 
the  high  explosive  in  the  flashing  process,  (2)  suggestions  of  explosive 
compositions  for  diagnostic  tests  of  flashing  or  other  degradation  mecha¬ 
nisms,  and  (3)  suggestions  for  practical  explosive  compositions  that  will 
minimise  flashing. 

•e  will  review  the  methods  of  calculating  detonation  parameters  and 
we  will  demonstrate  why  the  so-called  RUBY  code*  method,  although  by  no 
mesas  free  of  shortcomings,  appears  to  be  the  best  calculational  method 
currently  available.  Se  will  then  present  the  results  of  RUBY  calcula¬ 
tions  of  the  explosives  now  used  in  chemical  muni t ion*  and  calculations 
for  explosives  which  ms  believe  will  minimise  flashing.  As  this  pre¬ 
supposes  an  uaders*andlng  of  the  flashing  process,  some  justification 
of  the  hypotheses  adopted  for  the  flashing  process  will  be  appropriate 
at  that  stage  of  the  report.  Explosives  for  diagnostic  tests  of  the 
various  flashing  mechanisms  will  U>  proposed,  finally,  me  will  suggest 
a  practical  explosive  composition  that  should  mini alee  flashing  and  we 
will  present  some  recommendations  for  improving  detonation  calculations. 
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B.  General  Review  of  Detonation  Calculations 

Anyone  Involved  In  the  design  or  procurement  of  munitions  will,  at 
■>omc  tire,  be  faced  with  a  set  of  calculated  detonation  parameters.  Un¬ 
fortunately,  there  is  no  general  agreement  on  how  these  calculations 
should  be  made,  nor  is  there  general  agreement  in  the  results  of  the 
calculations.  Thus,  a  munition  designer  must  decide  for  himself  which 
calculations  he  should  believe.  Because  of  its  direct  bearing  on  our 
immediate  objective  of  preventing  agent  flashing,  we  will  present  here 
a  brief  state-of-the-art  review  of  detonatien  calculations. 

All  detonation  calculations  assume  equilibrium  conditions;  i.e. ,  all 
the  chemical  and  hydrodynamic  effects  are  assumed  to  have  reached  a  steady 
state.  Usually  calculations  are  carried  out  for  one-dimensional  detona¬ 
tions  so  that  the  complicated  effects  of  divergent  flow  can  be  neglected. 
In  common  practice  truly  one-dimensional  detonation  is  approached  only 
in  the  central  regions  of  a  sizable  explosive  charge,  at  a  considerable 
distance  from  the  point  of  initiation.  For  a  steady  one-dimensional  wave 
it  is  then  possible,  in  principle,  to  calculate  all  the  detonation  param¬ 
eters  (detonation  velocity,  particle  velocity,  pressure,  temperature, 
specific  volume,  and  product  composition)  from  conservation  of  mass,  con¬ 
servation  of  momentum,  conservation  of  energy,  an  equation  of  state  for 
the  detonation  products,  the  thermochemistry  of  products  and  reactants, 
and  the  Chapman- Jouguet  condition.  The  Chapman- Jouguet  condition  requires 
that,  at  equilibrium,  D  *  C  +  u,  where  D  is  the  detonation  velocity,  C  is 
the  local  sound  velocity,  and  u  is  the  particle  velocity.  Assuming  deto¬ 
nation  products  are  perfect  gases  and  neglecting  the  difference  between 
the  initial  pressure  and  tbe  final  pi^sure,  these  conditions  lead  to 


0  /o 
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V  +  1 
1 


-  1] 
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where  p  is  dmsity,  y  is  the  ratio  of  the  specific  heat  at  constant 
pressure  ana  constant  volume,  P  is  pressure,  Q  is  heat  of  reaction  per 
gram  of  ex-plosive,  T  is  temperature,  and  c  is  specific  heat  at  constant 
volume.  The  subscript  1  refers  to  the  equilibrium  detonation  state,  and 
the  subscript  0  refers  to  the  original,  unreacted  explosive. 

Other  equations  for  detonation  velocity  or  particle  velocity  are 
readily  obtainable,  but  for  the  moment  we  will  devote  our  attention  to 
the  three  equations  above,  since  they  serve  to  illustrate  how  a  detona¬ 
tion  calculation  proceeds,  how  the  calculation  is  improved  by  correcting 
for  gas  imperfections  (detonation  products  of  condensed  explosives  are 
anything  but  perfect  gases),  and  how  the  equations  can  be  used  for  making 
estimates  of  variations  in  the  detonation  parameters  produced  by  changes 
in  the  chemistry  of  the  process  (e.g.,  changes  in  the  heat  of  reaction  or 
the  number  of  moles  of  gas  in  the  products).  As  an  example,  consider  the 
detonation  of  acetylene  in  oxygen,  where  all  products  are  assumed  to  behave 
as  perfect  gases  but  there  is  insufficient  oxygen  to  produce  solely  Ha0  and 
C02 .  We  might  expect  the  net  reaction  to  be  something  as  follows: 

Ca H2  +  20,  -  COs  +  CO  +  H20 

If  this  equation  gives  the  actual  product  composition,  then  the  problem 
of  obtaining  detonation  parameters  is  nearly  solved.  The  known  products 
and  reactants  determine  ,  according  to 

Q  =2  (n  AH°)products  -  2  (n  AH°)reactants, 

1  i  1  f  j  J  1 

where  n  is  the  number  of  moles  per  gram  of  the  i-th  and  ,1-th  species  and 

AH°  is  the  standard  heat  of  formation.  The  known  products  determine  c, 
f 

and  Vj •  Thus,  all  the  necessary  information  is  available  to  solve  Eqs.  (l) 
through  (3).  Even  for  ideal  gas  behavior,  however,  the  products  are  not 
necessarily  those  given  by  a  balanced  chemical  equation.  For  our  example, 
we  would  certainly  have  to  consider  that  the  products  rearrange  according 
to  the  water  gas  equilibrium, 

C°s  +  Ha  =  Ha0  +  CO  (a) 
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which  is  temperature-controlled.  The  products  may  also  dissociate 
according  to 

C02  =  CO  +  1/2  Qa  (b) 

which  is  both  temperature-  and  pressure-dependent.  Thus  reaction  (a) 
leads  to 


p(h2o)p(co) 

kW  =  p(C02)p(H2) 


and  reaction  (b)  leads  to 


K(T.p) 


p(C0)p(02)1/2 

p(co2)pK2 


where  K  is  a  thermodynamic  equilibrium  constant,  p  is  partial  pressure, 
and  Pt  is  the  total  pressure  at  equilibrium.  In  making  the  calculation 
we  assume  an  equilibrium  temperature  and  pressure,  and  with  the  appropriate 
equilibrium  constants  for  these  assumed  conditions  we  calculate  the  pro¬ 
duct  composition.  Based  on  this  composition  we  obtain  Q1(  C x ,  and  and 
solve  Eqs.  (l)  through  (3).  Using  the  calculated  Tx  and  Px,  we  repeat 
the  equilibrium  calculations  and  use  the  new  composition  to  recalculate 
Tx,  Pj,  etc.  This  is  repeated  until  the  assumed  and  calculated  tempera¬ 
tures  and  pressures  are  brought  into  agreement.  It  should  be  realized 
immediately  that,  aside  from  any  other  effects  on  Eqs.  (l)  through  (3), 
gas  imperfections  will  affect  the  equilibrium  calculations  since  the 
partial  pressures  in  the  equilibrium  expressions  will  have  to  be  replaced 
by  fugacities2  (probably  unknown  for  the  regions  of  interest  of  products 
from  condensed  explosives).  In  addition,  the  total  pressure  terms,  such 
as  the  one  in  the  dissociation  of  002 ,  will  be  much  larger  for  the  highly 
imperfect  gases  than  for  ideal  gases. 

The  effect  of  gas  imperfections  can  be  even  more  drastic  on  Eqs.  (l) 
through  (3)  than  it  is  on  the  equilibrium  calculations.  One  of  the 
simplest  imperfection  corrections  that  can  be  made  is  to  assume  that  the 
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product  gases  are  rigid  spheres  occupying  a  definite  volume,  the  so- 
called  co-volume  a,  which  is  Independent  of  temperature  and  pressure. 

This  leads  to  the  Abel  equation  of  state  (a  simplification  of  Van  der  Waals' 
equation)  and  transforms  Eqs.  (l)  through  (3)  into 


v(  k,  +  o/v 
*  _  * _ o 

V  ~  k,  +  1 

o  i 


2Qt(kt  -  1) 


i  v  -  or 

o 


(4) 


(5) 


2Q,  kj 

1  ~  (kx  +  l)ct 


(6) 


where  v  =  1/p  and  kt  =  1  +  (ntR/cj)  =  Yx  if  0/  =  constant. 

Comparison  of  this  set  of  equations  with  Eqs.  (l)  through  (3)  shows 
that  if  product  composition  is  essentially  unaffected  by  the  gas  imper¬ 
fections  considered  (i.e„,  Qt  and  cx  are  essentially  the  same  for  Abel 
and  ideal  products),  then  the  only  parameter  that  is  appreciably  altered 

is  P.  in  the  sense  that  (P.)A.  ,  >  (P.)AJ  Theoretical  considerations, 

1  v  wAbel  '  17  ideal 

as  well  as  lack  of  agreement  between  calculated  and  observed  detonation 
parameters ,  make  untenable  the  assumption  of  an  a  independent  of  p  and  v. 
Several  attempts  have  been  made  to  improve  the  Abel  equation  by  removing 
the  restriction  that  a  be  constant.  Cook3  has  made  extensive  calcula¬ 
tions  using  a  PvT  equation  of  state  of  the  type 


Pv  =  nRT  +  o(v)P  (7) 

where  Q'(v)  is  obtained  from  experimental  measurements  or  detonation  ve¬ 
locity.  A  variation  of  the  form  of  Eq.  (7)  has  been  obtained  independently 
by  Paterson4  using  an  essentially  "first  principles"  (cell  theory)  approach 
rather  than  the  empirical  approach  of  Cook.  Paterson's  equation  of  state  is 


Pv  =  nRT[l  +  (b/v)  +  0.625(b/v)3  +  0.287(b/v)3  +  0.193(b/v)4]  (8) 
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where  b  =  £  n^,  and  bt ,  the  second  virial  coefficient  of  the  i-tfa  gas 
species,  is  assumed  constant.  Extensive  calculations  using  Paterson's 
method  are  griven  by  Taylor.5  Since  the  calculations  require  measured 
detonation  velocities  and  the  variation  of  detonation  velocity  with 
loading  density  of  the  explosive,  another  empirical  approach,  that  of 
Jones,®  is  based  on  an  equation  of  state  identical  with  Eq.  (7)  except 
that  or  =  Qf(p)  and  not  a  =  Qf(v). 

The  cell  theory  equation  of  state  of  Lennard-Jones  and  Devonshire, 
usually  called  LJD,  was  used  for  detonation  state  calculations  by  Fickett.7 
Ficket+.'s  treatment,  though  subject  to  some  criticism  on  theoretical 
grounds,  leads  to  certain  results,  confirmed  by  experiment,  that  are  not 
obtainable  by  other  calculational  methods  now  used. 

The  so-called  BKW  equation  of  state,  upon  which  RUBY  code  is  based, 
was  put  into  its  present  form  largely  by  Cowan  and  Fickett8  and  by  Mader.® 
It  Iji  semiempirical  in  nature  and  takes  the  form  of 


Pv 

RT 


1  +  Xe 


PX 


where 


X 


♦ z  °i°i 
v(T  +  6 )a 


(9) 


The  quantities  a,  (3 ,  and  6  are  adjustable  parameters  that  are  intended 
to  be  universal  constants,  t  is  a  scaling  factor  and  cr^  is  the  co-volume 
of  the  i-th  gas  species.  @  and  i  are  adjusted  to  make  the  calculated 
detonation  velocity  and  its  variation  m4th  packing  density  fit  the  ex¬ 
perimental  observation  on  a  "standard"  explosive.  Then  "a"  is  chosen 
to  reproduce  the  observed  pressure  of  the  "standard,"  Since  pressure 
would  go  to  infinity  when  T  -*  0,  a  factor  6  is  introduced  into  Eq.  (9). 
The  parameters,  a,  0,  t,  and  a,  have  been  frequently  reevaluated  and 
readjusted.  Thus  some  confusion  exists  among  calculations  done  at  dif¬ 
ferent  times  and  at  different  facilities.  Currently,  Mader9  uses  two 
"standards " — RDX  and  TNT — with  some  different  "constants"  for  each 
standard. 
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II  CALCULATIONS 


A.  Merita  and  Shortcomings  of  BKW  Equation  of  State 

Calculations  based  on  Abel's,  Cook's,  or  Paterson's  PvT  equations 
of  state  are  in  principle  applicable  only  to  thermodynamically  ideal  mix¬ 
tures  since  these  methods  neglect  the  energy  of  gas  imperfections.  It  is 
to  be  expected  that  detonation  products  of  condensed  explosives  ( "gases" 
such  as  CO2,  N2,  and  II2O  at  densities  greater  than  the  packing  density  of 
the  original  explosive)  are  far  from  thermodynamically  ideal.  As  these 
calculations  do  not  consider  any  repulsion  energy  of  the  detonation  prod¬ 
ucts,  the  detonation  temperatures  they  generate  will  be  too  high,  since 
the  entire  chemical  and  hydrodynamic  energy  increase  appears  as  an  in¬ 
crease  in  temperature  of  the  detonation  products.  Because  the  adjustable 
parameters  of  all  previously  mentioned  computational  procedures  were  chosen 
to  fit  experimentally  observed  detonation  velocities  (sometimes  also  det¬ 
onation  pressures),  and  in  view  of  the  fact  that  calculated  velocities, 
and  to  a  lesser  extent  pressures,  are  rather  Insensitive  to  the  form  of 
the  PvT  equation  of  state,  it  is  not  at  all  surprising  that  these  proce¬ 
dures  are  successful  in  predicting  detonation  velocities  and  detonation 
pressures.  A  much  more  severe  test  of  the  validity  of  a  computation  scheme 
is  to  be  found  in  comparing  calculated  detonation  temperatures  and  product 
compositions  with  experiment.  Unfortunately,  for  condensed  explosives 
there  are  few  (and  some  may  claim  no)  unambiguous  experimental  detona¬ 
tion  temperatures  and  certainly  no  product  compositions  measured  at  the 
detonation  equilibrium  state. 

Thus  the  choice  among  existing  PvT  equations  of  state  for  any  com¬ 
putational  scheme  should  be  made  primarily  on  the  basis  of  fundamental 
physics.  Consequently,  schemes  that  allow  for  intermolecular  energy 
should  be  considered  more  plausible  than  procedures  that  neglect  it. 
Fickett's  treatment  of  the  LJD  equation  of  state7  may  suffer  somewhat 
from  an  Inadequate  approximation  for  treating  gas  mixtures  (detonation 
products  are  usually  mixtures  and  not  single  species).  In  addition, 
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the  dependence  of  repulsion  energy  on  intermolecular  separation  may  not 
have  been  expressed  in  a  manner  that  agrees  best  with  theory;  nevertheless, 
his  treatment  qualitatively  predicts  the  shape  of  the  adiabatic  exponent 
(along  an  isentrope  through  the  equilibrium  detonation  state)  versus  ex¬ 
panded  volume  curve  actually  observed  for  several  explosives.10  Unfor¬ 
tunately,  Fickett's  method  is  not  retaily  adaptable  to  fairly  routine  cal¬ 
culations  . 

The  BKft'  equation  of  state  allows  only  forces  of  intermolecular 
repulsion  and  completely  neglects  any  forces  of  intermolecular  attrac¬ 
tion.  Thus,  the  BKW  equation  certainly  underestimates  detonatior  tem¬ 
perature  and  is  inadequate  for  isentrope  expansion  calculations;  however, 
in  its  RUBY  code  form1  it  does  lend  itself  to  routine  calculations.  The 
detonation  temperatures  and  the  product  composition  for  oxygen-poor  ex¬ 
plosives  obtained  from  RUBY  must  be  taken  with  a  considerable  amount  of 
skepticism.  It  is  possible,  however,  that  RUBY,  though  it  errr  in  absolute 
values  of  temperature  and  composition,  correctly  predicts  relative  values. 

In  any  event  RUBY  calculations,  with  all  their  shortcomings,  are  still  more 
plausible  than  computation?  based  on  the  several  variants  of  the  Abel 
equation. 

For  many  military  explosives  at  near  maximum  loading  density,  Mader9 
calculates  Tt  «  2500° K  while  Cook3  calculates  Tt  s*  5000°k  .  The  Pj's 
calculated  by  both  methods  are  roughly  the  same.  Let  us  examine  how  the 
differences  in  Tl  affect  equilibrium  product  composition.  Because  of  he 
high  detonation  pressures,  dissociation  reactions  of  the  type  C0a  =  C0+l/2  02 
will  be  repressed,  although  'ook  (because  of  his  higher  T*)  will  favor  dis¬ 
sociation.  In  the  so-called  soot  reactions, 


2C0 


CO,  +  c,  . 

(s' 


(c) 


CO  +  h2 


Ha0 


'(■) 


(«0 


high  pressure  favors  the  formation  of  solid  carbon.  The  following  tabula¬ 
tion  shows  that  for  a  given  pressure  the  soot  reactions  will  be  greatly 
enhanced  by  low  T1's: 


16 


Reaction 

Thermodynamic  Equilibrium  Constants 

At  2300 °K  At  5000°K 

(c) 

3.4  x  10-* 

1.3  x  10~T 

(d) 

2.7  x  10_a 

1.6  x  10~* 

Large  amounts  of  soot  are  commonly  observed  In  the  products  of  det¬ 
onation  of  military  explosives.  There  is  some  evidence  that  the  soot  is 
largely  in  the  amorphous  rather  than  in  the  graphitic  carbon  form.  It  is 
interesting  to  estimate  to  what  extent  calculated  detonation  product  com¬ 
positions  depend  on  the  form  of  the  solid  carbon.  Even  more  important 
would  be  estimates  of  how  incomplete  precipitation  of  solid  carbon  affects 
product  composition.  Such  Incomplete  precipitation  may  be  one  of  the 
most  common  causes  of  the  failure  of  an  actual  detonation  to  reach  thermo¬ 
dynamic  equilibrium.  Figure  1  presents  RUBY  calculations  for  Comp  B  under 
different  assumptions  for  the  fate  of  product  carbon.  It  is  obvious  from 
Fig.  1  that  the  form  of  the  solid  carbon — amorphous  or  graphitic — has 
little  effect  on  product  composition.  For  the  extreme  case  of  no 


U  ««*«< 


FIG.  1  PRODUCT  COMPOSITIONS  FOR  COMP  B-3  CALCULATED 
UNDER  THREE  ASSUMPTIONS  CONCERNING  CARBON 
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carbon  precipitation  the  products  become  richer  in  CH4  and  COj  primarily 
at  the  expense  of  H,0.  For  conditions  favoring  soot  reactions  (oxygen- 
deficient  explosives  at  high  packing  density) ,  one  can  make  fairly  ac¬ 
curate  estimates  of  the  detonation  product  composition  merely  by  allow¬ 
ing  all  the  hydrogen  atoms  of  the  original  explosive  to  form  Ha0,  with 
any  remaining  oxygen  combining  with  the  initially  present  carbon  to  form 
CO, .  Any  carbon  left  over  will  appear  as  solid  carbon  in  the  products 
and  any  nitrogen  atoms  originally  present  will  form  Na.  As  an  example, 
consider  RDX .  At  high  loading  density  the  reaction 

C3H€0#Ns  -  1.5  C0a  +  3  Ha0  +  3  N2  +  1.5 

is  much  favc-ed  over  the  reaction 

CjH.OgH,  -  3  CO  +  3  H,0  +  3  K, 

B.  Comparison  of  RUBY  with  Experiment 

It  has  already  been  mentioned  that  is  much  more  sensitive  to  the 
choice  of  a  PvT  equation  of  state  than  Pa .  Unfortunately,  completely 
unambiguous  experiments*  T, 's  are  not  available.  Experimental  tempera¬ 
tures  for  granular  explosives  must  be  considered  much  less  reliable  than 
experimental  temperatures  for  clear  liquid  explosives  because  of  "detona¬ 
tion  light"11  and  sir  shock  due  to  air  inclusions.  In  Tsble  I  published 

a 

experimental  valuea  of  Tj  and  P,  are  compared  with  Rader's  BKW  calcula¬ 
tions.  ( Kader  has  done  more  "curve-fitting"  and  has  published  calcula¬ 
tion  for  w  wider  variety  of  explosive  compositions  than  have  RUBY  users.) 
There  is  some  disagreement  In  experimental  temperatures  and  somewhat  less 
disagreement  in  experimental  pressures.  As  expected,,  measured  and 
calculated  pressures  agree  wtsll.  For  liquid  explosives,  the  temperatures, 
calculated  and  experimental,  agree  within  150°K  in  three  instances;  the 
lowest  experimental  values  are  at  least  250°K  higher  than  the  calculated 
la  four  cases,  and  in  two  cases  they  are  appreciably  lower  than  calculated 
values.  Urns,  when  allowance  is  made  for  the  uncertainty  in  measured  det¬ 
onation  temperature,  RUBY  calculations  and  experiments  are  in  reasonable 
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Table  X 

COMPARISON  OF  RUBY  CALCULATIONS  WITH  EXPERIMENTS 
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Table  I  ( Concluded) 


Note:  all  results  without  footnotes  as  well  as  all  calculations  (unless 
indicated)  from  llader.9 


1600  H 

(a)  Oxygen  balance  =  -  - - - - —  [2C  atoms  *■  —  atom  -  O  atom.] 

v  formula  wt  v  2 

(L)  A.  N.  Dremin  and  P.  F.  Pokhil,  Doklady  128(5),  989  (1959). 

(c)  I.  11,  Voskoboyr.lkov  and  A,  Ya.  Apits,  Doklady  130(4),  804  (1960). 

(d)  F.  C.  Gibson  et  a1.,  JAP  29,  628  (1958). 

(e)  J.  T.  A.  Burton  and  J.  A.  Hicks,  Nature  202 ,  75f  (1964). 

(f)  W.  C.  Davis,  B.  G.  Craig,  and  J.  3.  Ramsey,  "Failure  of  the  Chapaum- 

Jouguet  Theory  for  Liquid  and  Solid  Explosives,"  Preprints,  Fourth 
Symposium  ca  Detonation,  Vol .  I,  1965,  p.  A-50. 

( g)  D.  Price  and  H.  Hurwltz,  Ref.  12. 


agreement,  although  one  may  detect  a  tendenc;  for  RUBY  temperatures  to  be 
lower  than  observed  temperatures. 


C.  Comparison  of  Calculations  at  Different  Facilities 

A  a  already  mentioned,  the  -parameters  of  the  BSW  equation  •(  Eq .  9)  have 
undergone  some  charge  through  the  years.  Typical  variations  are  shown  in 
Tabic  II.  These  variations  result  in  the  differences  shown  in  Table  TII 
in  computed  detonation  properties  between  Los  Alr.jaos  Scientific  Laboratory- 
Naval  Ordnance  Laboratory  ( LASL-NOL)  and  Stanford  Research  Institute — 
Lawrence  Radiation  Laboratory  (SSl-UtL).  Generally,  the  ISL-SStl  T.  are 
lower  and  v,  are  higher  than  the  corresponding  LASL-ieOL  value*.  For 


T»ble  II 


BKV  PARAMETERS 


Parameter 

Mader* 

RDX  Fit 

TNT  Fit 

SSI  and  LRL 

a 

1/2 

1/2 

1/2 

s 

0.160 

0.0959 

0.100 

8  (°K) 

400 

400 

400 

1 

o-x  ( cc/kg) 

10.91 

12.69 

11.85 

H*0 

250 

250 

360 

CO, 

600 

600 

670 

CO 

390 

390 

; 

N, 

380 

380 

380 

H, 

180 

180 

180 

CH4 

528 

528 

528 

0, 

350 

350 

350 

_ 

comparison  this  table  also  includes  Cook's  (based  on  Eq .  7),  Paterson’s 
(based  on  Eq .  8),  and  Fickett's  results  (based  on  tbe  LJD  equation  of 
state).  Compared  to  tbe  Under  ( LASL)  or  SOL1*  calculations  Oook  and 
Paterson  obtain  sosewhat  lover  pressures  and  such  higher  temperatures. 
Fickett's  results  for  Px  are  also  lower  than  Msder’s  and  hi*  Tx's  are 
between  those  of  Mader  and  Cook.  Equilibrium  product  compositions  ob¬ 
tained  by  Under  are  compared  to  our  calculations  in  Table  IV.  Generally 
Bader's  results  show  a  little  more  solid  carbon  than  ours  (  at  .<»  expense 
of  CO.  and  CC  and  a®re  H,0  (at  the  expense  of  minor  products  such  as 
CH4  and  RBj  ) . 


Table  III  (Concluded) 


Calculated  by 

Explosive 

B 

D 

( mm/M-sec) 

JSBj 

T i 

(°K) 

vi 

(cais/g) 

LASL  (BKW  Code) 

RDX 

1.80  j 

8.754 

347 

2590 

0.416 

NOL  RUBY 

8.567 

341 

2670 

0.412 

LRL 

- 

- 

- 

- 

SRI  f 

8.791 

353 

0.523 

Cook 

8.84 

?7S 

0.450 

LJD-Fickett 

8.798 

324 

- 

Paterson 

1 

1 

f 

- 

- 

- 

- 

LASL  (BKW  Code) 

1  | 

PETN 

1.77 

8.421 

319 

2830 

0.421 

NOL  RUBY 

8.337 

324 

2810 

0.414 

LRL  j 

1 

- 

- 

- 

- 

SRI  f 

8.575 

2690 

0.413 

Cook 

8.47 

260  1 

i 

5630 

0.453 

LJD-Fickett 

- 

- 

- 

- 

Paterson 

r 

\ 

r 

I 

- 

- 

- 

- 

LASL  (BKW  Code) 

PETN 

1.20 

6.599 

3590 

0.603 

NOL  RUBY 

1  ; 

i 

- 

LRL  1 

i 

I 

- 

- 

SRI  ▼ 

6.747 

3440 

0.552 

Cook 

| 

6.40 

5310 

0.630 

LJD-Fickett 

- 

- 

- 

Paterson 

r 

r 

6.60 

118 

5200 

0.620 

LASL  (BKW  Code) 

HMX 

1.90 

9.159 

395 

2360 

0.396 

NOL  RUBY 

8.895 

382 

2510 

0.392 

LRL  | 

9.084 

388 

2140 

0.374 

SRI  f 

£.110 

390 

2270 

0.373 

Cook 

- 

- 

- 

l 

LJD-Fickett 


Paterson 
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Table  IV 


CALCULATED  PRODUCT  COMPOSITIONS 
(mole  $) 


Product 

HMX 

(p0  =  1.6  g/cc) 

...  , , 

TNT 

(p0=1.6  g/cc) 

TOT 

(p0=1.2  g/cc) 

RDX 

(p0»1.8  g/cc) 

LASL 

SRI 

LASL 

SRI 

LASL 

SRI 

LASL 

SRI 

lo  ,9 

17.8 

P  1 

H.4 

17.1 

16.5 

18.0 

mmmMm 

M  M 

1.5 

2.1 

■  :  1; 

9.5 

7.7 

0.2 

1.0 

ch4 

- 

m 

- 

. 

- 

3.2 

- 

0.8 

C,  \ 

(a) 

15.9 

14.2 

46,7 

44.1 

42.8 

40.2 

16.5 

14.4 

HaO 

33.3 

30.6 

22.7 

17.7 

22,4 

16.6 

33.3 

31.3 

h2 

- 

- 

- 

- 

0.2 

0.3 

- 

- 

nh3 

- 

1.2 

- 

0.7 

- 

0.7 

- 

0.9 

n2 

33,3 

32.2 

13.6 

14.1 

13.6 

14.3 

33.3 

33.7 

D.  Calculations  for  Comp  B  and  Comp  B-3/Oxamide 

The  dissemination  charge  presently  used  is  Comp  B-3.  Consequently, 
it  is  important  to  estimate  its  detonation  parameters  and  product  composi¬ 
tion.  Using  the  input  parameters  of  Table  V,  the  results  of  the  calcula¬ 
tion  are  shown  as  functions  of  loading  density  in  Figs.  2  and  3. 

It  has  been  suggested1 3  that  an  explosive  charge  consisting  of  90/10 
weight  percent  (all  subsequent  compositions  will  be  given  in  weight  per¬ 
cent)  of  Comp  B/cxamide  could  minimize  flashing  because  is  expected 

to  be  less  for  this  mixture  than  for  Comp  B.  That  this  appears  to  be 
so  is  also  shown  in  Figs.  2  and  3.  Note,  however,  that  the  product 
composition  of  Comp  B  is  .he  same  as  the  product  composition  of  90/10 
Comp  B/oxamide.  Thus  the  combustibility  of  the  products  for  the  two 
explosives  is  the  same,  but  the  lower  temperature  can  be  expected  to 
reduce  the  probability  of  product  ignition  in  the  case  of  the  oxamide 
mixture.  The  90/10  Comp  B/oxamide  composition  does  not  eliminate 
flashing.14 
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LOADING  DCNSITY - 9/ao5 

FIG.  2  CALCULATED  DETONATION  PARAMETERS  FOR  COMP  B'3 
AND  90/10  COMP  B-3/DXAMIDE-  TK#  solid  lmo«  or*  for  pur# 
Com?  B'3  and  H»o  dasfvod  linos  aro  far  Comp  B-J  axoietdo. 


.8  1.0 


I  '  1.2  1.4  1.6 

LOADING  DENSITY - fl/cm3 


I.Z  l.e  i.o 

LOADING  DENSITY -  fl/cm3 

T»-«»oo-ro 


FIG.  3  DETONATION  PRODUCTS  OF  COMP  B~3  AND  90/10  COMP  B-3/0XAMIDE 


£.  Calculations  for  TWT  and  Amatols 

Flashing  could  be  caused  by  the  ignition  of  an  agent -air  Mixture  by 
burning  detonation  products.  If  the  detonation  products  are  incombustible 
this  cause  of  flashing  is  obviously  eliminated.  An  explosive  formulated 
to  produce  solely  CO,,  H,0,  and  N,  is  said  to  be  oxygen-balanced  and 
produces  no  combustible  products. 

Most  military  explosives  have  negative  oxygen  balance;  i.e.,  they 
contain  insufficient  oxygen  to  convert  all  their  carbon  atoms  to  CO,  and 
all  their  hydrogen  atoms  to  steam.  In  industrial  explosives  a»onium 


nitrate  (AN)  is  commonly  ur«d  as  an  ene  -getic  oxygen  supplier.  Ammonium 
nitrate  has  several  drawbacks  such  a?  hygrorcopici ty (  a  tendency  to  cake, 

and  a  phase  change  at  around  32°C  which  produces  a  large  volume  increase. 

'  -•»_  - - x.  .....  . 


Nevertheless,  certain  well-known  mixtures  containing  AN  could  be  valuable 
in  diagnostic  tests  of  flashing  mechanisms.  It  is  for  this  reason  that 
we  have  made  calculations  for  TNT,  50/50  TNT/AN,  and  20/80  TNT/ AN.  The 
TNT/ AN  mixtures  are  called  amatols. 


Since  TNT  is  very  deficient  in  oxygen,  a  large  proportion  of  its 
products  will  be  combustible  in  air.  The  50/50  amatol  has  slightly  nega¬ 
tive  oxygen  balance  while  the  80/20  amatol  is  very  slightly  oxygen-rich. 
The  results  of  the  calculations  (see  Tables  V  and  VI  for  input  data)  are 
shown  in  Figs.  4  and  5.  Rather  surprisingly,  Tj  for  the  amatols  is  cal¬ 
culated  to  be  substantially  lower  than  that  for  the  TNT.  In  part,  this 
is  a  consequence  of  the  great.ei  thermodynamic  stability  of  the  An  as  com¬ 
pared  to  TIT.  (more  negative  A  H°  ;  see  Table  VI),  which  results  in  a 
smaller  hett  of  reaction,  Q1(  for  the  amatols  than  for  TNT.  Possibly 
even  more  Important  is  the  fact  that  Zn1cl  increases  substantially  in 
going  from  TNT  to  amatol.  This  also  lowers  Tx ,  since  in  the  limit  of  a 
perfect  gas  (Eq.  3)  Tt  varies  as  Q1/En1c1.  The  substantial  Increase  in 
calculated  pressure  of  amatols  over  that  of  TNT  (compared  at  a  fixed  p0) 
is  puzzling.  The  interrelation  of  calculated  parameters,  e.g. ,  the  de¬ 
pendence  of  Pj  on  Qj,  will  be  examined  in  b  later  section. 


F.  Calculations  of  HMX  and  HMX/AN 

For  practical  reasons  as  well  as  for  the  further  diagnosis  of  flash¬ 
ing,  it  may  be  desirable  to  have  an  explosive  which  produces  a  much  higher 
detonation  pressure  than  the  amatols  but  like  the  20/80  amatol  contains 
no  combustible  products.  As  can  be  seen  in  Figs.  6  and  7  this  objective 
appears  attainable  with  a  52/48  HMX/AN  mix,  which  we  call  "Ametex."  Again, 
as  for  TNT  and  amatols,  Tx  for  HMX/AN  is  much  lower  and  Pt  is  slightly 
higher  than  the  corresponding  values  for  HMX.  This  \\nexpected  result  will 
be  discussed  later.  The  well-known  explosive  PSX  9404  consists  of  94* 

HMX  and  6*  binder. 
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DETONATION  PRODUCTS  OP  TNT,  50/50  AMATOL,  AND  80/20  AMATOL 


m 


FIG.  6  CALCULATED  DETONATION 

PARAMETERS  FOR  HMX  AMO  52  « 
HMXAN.  TV#  *olid  lin*i  «  fef  pgr* 
HMX  mmd  #*»  4c*te*<j  iin**  or*  (or  52.  48 
HMX  AN  - 
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FIG.  '  DETONATION  PRODUCTS  OF  HHX  AND  52  48  HMX  AN 


C ,  Appltcsbiil tv  of  i  tooatlon  State?  Calculations  to  Flashing 

At  this  point  it  sill  be  instructive  to  rwiiaine  the  product  com¬ 
position  versus  p6  corses  far  oxygen -ha  la  need  explosives  '  F  x  gs .  5  sad  7) 
st  is  obvious  that  product  composition  for  oxygen -balanced  explosives  is 
550t  influenced  by  *b@  original  loading  density  Pe,  of  the  explosive. 


Because  p0  strongly  influences  Pt ,  and  to  a  lesser  degree  Tj ,  one  Bust 
conclude  that  product  coBposition  for  oxygen -balanced  explosives  is  nearly 
Independent  of  Tj  and  Pj .  ITiis  Is  a  ver®  important  conclusion  since  it 
implies  that  product  compositions  calculated  for  the  detonation  state  of 
oxygen-balanced  explosives  will  remain  essentially  unchanged  as  these 
products  expand. 

Flashing  probably  occurs  during  or  after  detouation  product  expansion. 
Thus  there  is  always  the  question  of  whether  detonation  state  calculations 
are  applicable  to  the  flashing  problem.  This  question  is  answered  in  the 
affirmative  for  oxygen -balanced  explosives.  For  oxygen -poor  explosives  an 
unequivocal  answer  is  not  possible.  RUBY  product  composition  calculations 
along  an  expansion  isentrope  passing  through  the  detonation  equilibrium 
state  are  not  very  meaningful  since  the  BKW  equation  must  fail  as  the 
pressure  decreases  (see  Section  IIA).  Nevertheless,  such  calculations 
have  been  made  at  LASL  and  at  LRL,  ad  they  suggest,  as  might  be  expected 
from  applying  Le  Chatelier's  principle  to  the  chemical  equations  D  and  c. 
that  the  main  composition  change  dur  tig  expansion  (pressure  decrease)  is 
a  decrease  in  solid  carbon.  For  sufficiently  great  expansion  there  is  cf 
course  a  temperature  drop  ns  well  as  .1  pressure  decrease.  The  temperature 
decrease  irplies  that  at  seme  stage  of  the  expansion  ^he  various  chemit  :1 
equilibria  will  "freeze"  because  the  temperature  of  the  system  becomes  too 
low  to  maae  further  chemical  readjustments  sufficiently  rapid. 

Ornellas  et  al.1*  at  LRI.  have  estimated  the  temperature  rep  on  below 
which  equilibria  are  "frozen"  to  be  150C  lBOOck,  and  they  have  measured 
the  prod-c  composition  of  PETS  detonated  In  a  calorimeter  r  *  cooled  to 
25  °C.  In  one  set  of  experiments  the  PfcTN  was  confined  so  that  in  ‘jing 
work  on  the  confinerent  the  temperature  dropped  rapiolv.  In  an  unconfined 
snmple,  beet-use  of  reverberations  at  the  calorimeter  walls,  the  temperature 
dropped  much  nore  gradually  than  in  the  comined  sample.  Their  results 
are  coapa’-ed  with  ESV  calculations  in  Table  HI.*  It  is  apparent  that  the 
isentrope  calculation  '  for  "freeze-out"  teftperatures  between  1500  arwr 
1 80C )  does  net  agree  well  •  1  th  experiment.  In  fact  it  may  be  argued 

* 

*e  do  not  kco*  of  any  other  reliable  product  composition  mtarurveenti- 

other  than  the  ones  in  Xcf.  15. 
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COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL  EXPANDED 
PRODUCT  COMPOSITIONS  FOR  PETN 
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that  a  detonation  state  calculation  for  low  packing  density  PSTN,  also 
shown  in  Table  VIZ,  agrees  with  experiment  about  as  well  as  the  isentrope 
calculation.  All  this  still  leaves  unanswered  the  question  of  the  ap¬ 
plicability  of  detonation  state  calculations  (for  oxygen-negative  ex¬ 
plosives)  to  the  expanded  state.  It  suggests  that  a  crude  estimate  of 
expanded  product  composition  may  be  made  from  detonation  state  calcula¬ 
tions  for  explosives  at  low  packing  density, 

H.  Calculations  for  HMX  and  HMX/AP 

In  Section  IIP  it  was  shown  that  an  energetic  explosive  containing 
no  combustible  products  can  be  made  by  mixing  HMX  and  AH.  IM3  mixture, 
however,  would  be  unsuitable  for  military  use  because  of  the  poor  storage 
qualities  of  AN.  It  is  believed  that  a  practical  military  explosive  c  n 
be  made  from  mixtures  of  HMX  and  ammonium  perchlorate  (AP)  imbedded  in 
a  suitable  plastic  matrix.  Calculations  for  the  oxygen-balanced  61.2/38.6 
HMX/AP  mixtures  are  shown  in  Figs.  8  and  9.  Comparing  these  results  with 
the  calculations  for  HMX,  given  in  Figs.  6  and  7,  shows  that  Pj  is  greater 
and  Tj  is  less  for  HMX/AP. 

I .  Flashing  Mechanism 

Several  investigations14 *18 i1T  suggest  that  flashing  results  pri¬ 
marily  from  Lb'  Ignition  of  an  agent-air  mixture  by  burning  of  detona¬ 
tion  products.  Combustible  rixtures  in  air  for  either  products  or  agent 
are  formed  only  after  the  whole  system  has  undergone  a  considerable  degree 
of  expansion.  The  energy  transport  from  the  hot  burning  products  to  the 
still  unignited  agent-air  mix  is  presumably  via  conduction,  convection, 
and  radiation,  all  of  which  are  slow  compared  to  hydrodynamic  transport 
phenomena.  The  detonation  temperature  is  rapidly  quenched  if  detona¬ 
tion  products  expand  and  do  work  on  their  surroundings.  Thus  ignition 
of  agent-air  mixtures  by  heat  transport  is  most  unlikely  if  the  only 
source  of  high  temperature  is  ,  However,  afterburning  of  detonation 
products  in  air  can  maintain  regions  of  high  temperature  (of  the  order 
of  Tj  or  even  higher)  for  times  sufficiently  long  to  allow  appreciable 
heat  transport  to  the  agent-air  mix.  Detonation  products  that  are  in¬ 
combustible  can  no  longer  act  as  an  energy  stimulus  of  sufficiently  long 
duration  to  ignite  the  agent-air  mixture. 
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FIG.  8  CALCULATED  DETONATION  PARAMETERS  FOR 
61.2/38.8  HMX/AP 
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LOADING  DENSITY - g/C*5 

•  A.  *900  -  i«« 

FIG.  9  DETONATION  PRODUCTS  OF  61.2/38.8  HMX'AP 


There  appears  to  be  Insufficient  tine  for  enough  heat  transfer,  froai 
essentially  unexpanded  detonation  products  to  the  agent,  to  cause  flashing. 
In  any  event  nost  agents  per  se ,  and  not  as  agent-air  alxtures,  weald  be 
difficult  to  ignite.  However,  if  Tt  Is  very  large,  it  la  conceivable  that 
enough  heat  can  be  transferred  to  the  agent  to  cause  som  of  1*  to  degrade 
(e.g.,  pyrolixe),  without  eeittlng  light,  before  7,  Is  drastically  reduced 
by  expansion. 

There  is  yet  another  possible  agent  degradation  eechanise  in  which 


the  shock  sent  into  the  agent  by  tbe  uetooating  explosive  creates  hot 
regions  of  sufficiently  high  temperature  to  cause  degradation. 


It  may  be  desirable  to  assess  the  relative  importance  of  these 
degradation  schemes.  To  that  end  ve  have  listed  in  Table  VIII  various 
explosive  compositions  that  could  be  used  for  diagnostic  teats.  The 
choice  of  the  explosives  was  made  on  the  basis  of  our  calculations  as 
well  as  published  calculations  of  others.  All  explosive  parameters  are 
rated  relative  to  Comp  B,  the  presently  used  dissemination  charge. 

\ 

J .  Interrelation  Among  Calculated  Petonation  Parameters 

It  is  most  desirable  to  be  able  to  check  the  internal  consistency 
of  the  results  of  an  involved  computational  procedure  such  as  RUBY.  Such 
checks  should  be  capable  of  pointing  out  errors  in  input  data  and/or 
machine  malfunctioning.  The  most  obvious  and  least  ambiguous  check  is 
naturally  that  of  mass  balance;  i.e.,  all  atoms  originally  present  must 
appear  in  the  calculated  products.  A  correct  mass  balance  is  a  necessary 
but  not  a  sufficient  condition  for  the  internal  consistency  of  the  calcula¬ 
tion.  Other  checks  are  needed  in  the  usual  caae  where  there  are  no  experi¬ 
mental  data  to  compare  with  the  calculation. 

As  shown  in  the  appendix,  a  simple  relation  between  Qt  and  Pj  is 
obtained  if  the  EPv  equation  of  state  is  E  =  Pv^r  -  1,  where  T  is 

-(&tn  P^dtn  v)g.  This  relation  is: 

~  SPoQjfr  -  1)  (10) 

Many  calculations*’1*  give  T  ~  constant  ~  3.  Thus,  according  to  Eq.  (10), 
a  plot  of  Pi  versus  P0Qi  should  be  s  straight  line  of  slope  2(T  -  l). 
Experimental  and  calculated  P,'s  are  plotted  versus  computed  p0Q,'s  in 
Fig.  10.  Actually  the  slightly  sore  accurate  p0  [Qt  +  njRT0  (T/r  -  l)] 
was  used  rather  than  P0Qi<  Most  of  the  points  lie  along  a  line  of  the 
expected  slope,  2(f  -  1),  with  P  -  3,  but  the  line  does  not  pass  through 
the  origin.  This  implies  that  Eq,  (10)  should  contain  an  additional  con¬ 
stant  and  that  the  simple  EPv  equation  of  state  used  in  the  appendix  is  not 
consistent  with  the  BKW  equation  of  state. 

In  Section  1IA  it  was  shown  that,  if  P,  is  large,  product  com¬ 
position  and  consequently  Qt  can  be  estimated  fairly  accurately  strictly 
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COMPARISON  OF  POSSIBLE  DISSEMINATION  TEST  EXPLOSIVES 
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FIG.  10  C0RRE1  MION  BETWEEN  DETONATION  PRESSURE  AND  HEAT  OF  REACTION 
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on  the  basis  of  bass  balance  and  thermochemistry,  l.e.,  without  recourse 
to  any  *vr  equation  of  state.  Thus  the  plot  of  Fig.  10,  regardless  of 
its  theoretical  significance  or  lack  of  it,  provides  a  method  of  corre¬ 
lating  most  of  the  computed  and  experimental  Px  -  (^  results.  For  the 
LASL-NQL  computations,®  »i3  only  THM  and  HX  do  not  fit  the  proposed  corre¬ 
lation.  Both  of  these  materials  are  oxygen-rich  and  have  small  (J^’s. 

Other  ox/gen-rich  materials  such  as  BTNEN,  NG,  and  TC  (see  Fig.  10)  fall 
close  to  the  proposed  correlation.  In  our  calculations  all  oxygen-balanced 
or  slightly  oxygen-rich  materials  have  a  higher  ?1  than  the  proposed  corre¬ 
lation  would  predict  for  their  respective  Q^'s.  Since  we  expect  all  Qx’s 
to  be  reasonably  correct,  this  could  Indicate  that  our  RUBY  computations 
for  oxygen-balanced  mixtures  somehow  give  P1  ' s  that  are  too  high— a  con¬ 
clusion  that  is  not  entirely  unexpected.  Ifce  necessary  "curve-fitting" 
for  RUBY  and  BKW  codes  was  done  using  the  oxygen-poor  RDX  and  TNT  as 
"standard"  explosives.  Hue  equation  of  state  parameters  (f,  0,  etc.) 
from  this  fit  may  not  be  applicable  to  oxygen-balanced  mixtures.  Moreover, 
the  co-volumes  (or^  of  Eq.  9)  for  1^,0  and  CCl,  used  in  our  calculations  are 
larger  than  those  used  by  LAFL  and  NOL  (see  Table  II).  Consequently,  the 
differences  between  SRI  and  LASL-NQL  calculations  for  oxygen- balanced 
explosives,  whose  products  contain  large  proportions  of  ^0  and  CCfc  ,  may 
reflect  this  co-volume  difference. 
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Ill  CONCLUSIONS  AND  RECOMMENDATIONS 


If  flashing  is  caused  by  ignition  of  agent-air  Mixtures  by  after¬ 
burning  detonation  products  of  the  dissemination  charge,  then  flashing 
can  be  elininated  by: 

1.  Formulating  the  explosive  charge  to  produce  incombustible 
products. 

2.  Designing  the  munition  so  that  detonation  products  expand 
symmetrically.  IXiring  expansion  the  products  do  work  on 
the  surroundings  and  are  cooled  before  they  come  into  con¬ 
tact  with  the  air.  It  has  been  shown14  that  Comp  B  charges 
do  not  flash  if  they  are  synetrically  surrounded  by  agent 
simulant  and  symmetrically  detonated. 

3.  Using  explosives  of  low  detonation  temperature  to  help  re¬ 
duce  flashing.  This  approach  appears  to  be  much  less 
promising  than  1  or  2.  In  any  event  it  appears  desirable 
to  keep  detonation  temperatures  low.  Even  if  no  flashing 
is  observed  a  high  detonation  temperature  could  cause 
localized  degradation  of  the  agent. 

The  easiest  and  surest  way  to  formulate  explosives  producing  incom¬ 
bustible  products  is  to  include  sufficient  oxygen-containing  compounds 
in  the  original  explosive  to  convert  all  the  carbon  atoms  and  hydrogen 
atoms  of  the  explosive  to  CC^  andH,0.  Excess  oxygen  should  be  avoided 
as  it  may  interact  with  the  agent.  Ideally  the  oxygen-balanced  explosive 
should  be  a  single  compound,  e.g.,  ethylene  glycol  dinitrate,  since  this 
avoids  all  problems  of  product  nonuniformity  due  to:  (1)  nonuniform 
and/or  insufficiently  intimate  mixing  of  explosive  constituents  and  (2) 
different  chemical  reaction  rates  for  the  various  compounds  in  the  mix¬ 
ture.  Unfortunately,  for  military  uses  the  requirement  of  an  oxygen- 
balanced  single  compound  is  too  restrictive.  It  seems  better  to  con¬ 
sider  oxygsn-bslanced  mixtures  but  to  exercise  great  care  in  selecting 
and  preparing  them. 
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The  preceding  sections  have  demonstrated  that  the  products  of  an 
oxygen- balanced  explosive,  provided  equilibrium  is  achieved,  are  not 
affected  by  Initial  density,  expansion  of  the  products,  or  even  the 
method  of  calculation.  The  other  computed  detonation  parameters,  such 
as  pressure  or  temperature,  of  course  depend  on  the  method  of  calculation. 
At  present,  any  computed  detonation  temperatures  must  be  accepted  with 
considerable  skepticism.  However,  viewed  in  terms  of  the  fundamental 
physics  involved,  temperature  calculations  based  on  the  BKW  equation  of 
state  are  inherently  more  plausible  than  calculations  based  on  modified 
Abel  equations  of  state.  The  present  RUBY  code  is  too  inflexible  to 
permit  the  use  of  other  equations  of  state  than  BKW,  which,  as  has  been 
discussed,  has  serious  theoretical  shortcomings.  A  very  general  computa¬ 
tional  code,  now  being  completed  at  SRI,  will  be  much  more  flexible  than 
RUBY.  This  new  code  should  be  employed  to  develop  a  theoretically  sound 
eq  ntion  of  state,  to  test  ic  against  available  data,  and  to  suggest 
areas  "'ssibly  temperature  measurements)  where  new  data  must  be  obtained. 

Shock  heating  of  the  agent  by  the  detonating  dissemination  charge 
may  cause  localized  agent  degradation.  RUBY  usually  provides  reliable 
estimates  of  detonation  pressure.  Thus  from  RUBY  results  one  can  suggest 
explosives  for  diagnostic  tests  of  the  importance  of  shock  heating 
(Table  VIl). 
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GLOSSARY 


b  Second  virial  coefficient 

C  Local  sound  veloc'ty 

c  Specific  heat  at  constant  volume 

D  Detonation  velocity 

E  internal  energy 

o 

AKf  Standard  heat  of  formation 

X  Thermodynamic  equilibrium  constant 

n  dumber  of  moles  per  gram  of  explosive 
P  Total  pressure 

p  Partia  pressure 

Q  Heat  of  reaction  per  gram  of  explosive 

R  Universal  gas  constant 

T  Temperature,  °K 

u  Particle  velocity 

v  Specific  ^lume 

a  Co- volume 

Y  Ratio  of  specific  hea‘  constant  pressure  ana  constant  volume 

p  Density 

V  bertiing  factor 

Subscripts 

,'}  Original ,  unreacted  explosive 

1  Equilibrium  detonation  state 

S  Constant  entropy 
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APPENDIX 


From  the  conservation  equations  for  a  one-dimensional  shock, 

Ei  '  E°  =  I  Vv<>  "  V  (A"1) 

if  Pj  »  P0.  The  subscript  j  refers  to  the  equilibrium  detonation  state 
(the  Chapman-Jouguet  state),  the  subscript  0  refers  to  the  unreacted  ex¬ 
plosive  at  initial  temperature  and  pressure,  and  the  subscript  1  refers 
to  equilibrium  detonation  products  which  are  not  necessarily  at  P  and 

J 

T^.  For  explosive  (T0,F0)  -*  equilibrium  detonation  products  (T0,P0), 

AE  =  Qi  +  niRT0  =  Eq ( T0 , P0 )  -  E^To.Po)  (A-2) 

Since  Ej  =  E^T  ,P  )  and  Ej  =  E0(T0,Po),  Eq.  (A-2)  gives 

Ej  -  E0  =  Bt(T  P  )  -  E,(T0,Po)  - 

and  substituting  in  Eq.  (A-l)  results  in 

I  Pj(vo  '  Vj)  =  Ei(Tj-Pj)  -  Ei(To,Po)  -  AE  (A-3) 

The  Chapman-Jouguet  condition  is  expressed  by 

v^/v0  =  Po/Pj  =  r/(r  +  1)  (a-4  ) 


where 


r 


/ 

A 


£n  p\ 

s 


•19 


For  the  so-called  polytropic  gas,  at  constant  entropy  S,  P  *  Ap1 
where  A  ■  A(S).  Under  certain  conditions1* > 18  this  leads  to  the  fol¬ 
lowing  EPv  equation  of  state 


Combining  Eqs.  (A-3l,  (A-4),  and  (A-5)  results  in 


(A-S) 


(A-6) 


If  the  detonation  products  at  (P0  .TJ,  )  are  perfect  gases,  Eq.  (A-S)  is 
further  simplified  to 


Ae  + 


nx  RT 
V  -  1 


Qi  +  nx  RTi  (1  + 


Qi 


n,  RT0  V 

V^nr 


since  >:>  RT0y/y  -  1  in  most  cases.  Thus  for  a  consistent  set  of 
dimensions  (cirs)  a  plot  of  P  versus  p0Q1  should  be  a  straight  line  of 
slope  2(F  -  1),  sin«.e  I'  -  constant  ~  3. 
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